The decay of excited He from Stochastic Density-Functional Theory: a quantum 

measurement theory interpretation 
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Recently, time-dependent current-density functional theory has been extended to include the 
dynamical interaction of quantum systems with external environments [Phys. Rev. Lett. 98, 
226403 (2007)]. Here we show that such a theory allows us to study a fundamentally important 
class of phenomena previously inaccessible by standard density-functional methods: the decay of 
excited systems. As an example we study the decay of an ensemble of excited He atoms, and discuss 
these results in the context of quantum measurement theory. 



Density- functional theory (DFT), in both its ground- 
state and time-dependent versions [l|, l^j y, Ij, [Sj has 
become the method of choice to study several equilib- 
rium and non-equilibrium properties of interacting many- 
particle systems evolving under Hamiltonian dynamics. 

There is, however, a large class of physical problems 
where the dynamical interaction of a quantum system 
with an external environment needs to be taken into ac- 
count. To this class of open quantum systems belongs also 
one of the most basic tenets of Quantum Theory, namely 
the non-unitary evolution of a quantum state due to the 
measurement by an apparatus. Non-unitary quantum 
evolution also pertains to processes where the energy of 
the quantum system relaxes into the degrees of freedom of 
a bath or reservoir, like, e.g., the decay of exited systems. 
An understanding of such processes from a microscopic 
point of view would represent a substantial advancement 
in the study of open quantum systems. 

To address the above issues, Di Ventra and D'Agosta 
(DD) have recently proved [6| that given an initial quan- 
tum state, and an operator V that describes the inter- 
action of a many-body system with an external environ- 
ment, two external vector potentials A(r,t) and A'(r,t) 
that produce the same ensemble-averaged current den- 
sity, j(r, i), must necessarily coincide, up to a gauge 
transformation. The DD theorem thus extends the pre- 
vious theorems of dynamical DFT (that are one of its 
corollaries ii V = 0), and allows for the first- principles 
description of the dynamics of open quantum systems 
via effective single-particle equations. This theory has 
been named Stochastic Time-Dependent Current-DFT 
(Stochastic TD-CDFT). 

Here we apply the above theory to a previously inac- 
cessible problem via standard DFT methods: the decay 
of an ensemble of excited He atoms. In addition, we inter- 
pret the results in the context of quantum measurement 
theory by showing that the interaction with the environ- 
ment can be viewed as a continuous "measurement" of 
the state of the system, thus making a connection be- 
tween density-functional theory and quantum measure- 
ment theory. 



We consider two cases: 1) an ensemble of excited He"*" 
atoms, whose dynamics can be directly compared with 
the one obtained from a density-matrix approach. 2) An 
ensemble of neutral excited He atoms. Our results reveal 
unexpected features of this problem, like the dampening 
and modification of high-frequency oscillations during en- 
ergy relaxation of the ensemble towards its ground state. 

The starting point of Stochastic TD-CDFT is the 
stochastic equation of motion of an auxiliary Kohn-Sham 
(KS) Slater determinant v]/^-5 built out of single-particle 
KS states (j)a (atomic units are used throughout this pa- 
per) 
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with A{fi,t), an arbitrary external vector potential, 
Axc[j{j^,t),\'io),V] the exchange-correlation vector po- 
tential (which is a functional of the average current j, the 
initial condition |^o)i and the operator V), and VniJ^jt) 
the Hartree potential. The quantity r has dimensions 
of time. Without loss of generality the stochastic pro- 
cess, £{t), is chosen such that it has both zero ensemble 
average and J— autocorrelation, i.e. [20|, 



e(t) = 0; e{t)i{t') = TS{t - 1') , 



(3) 



where the symbol '^^ indicates the average over a statisti- 
cal ensemble of identical systems all prepared in the same 
initial quantum state |\l/o)- For the particular choice of 
bath operator we will make in this paper (Eq. ^), which 
acts on single-particle states only, the stochastic equa- 
tion ^ is simply 

dtM-t) - -iHKsMt) - ^v^vMt) + eit)vMt) , (4) 

where a contains also the spin degrees of freedom. 



The use of a stochastic Schrodinger equation in the 
context of DFT, and not of an equation of motion for the 
density matrix, is because in DFT the KS Hamiltonians 
depend on the density (and/or the current density), and 
therefore they are, in general, different for the different 
elements of the ensemble. This does not generally guar- 
antee a closed equation of motion for the single-particle 
KS density matrix of the mixed state [0|. 

As mentioned previously, our aim is to describe the 
decay of excited electrons bound to a He nucleus. The 
electrons are prepared in some initial excited state, and 
evolve into the ground state as a result of the stochastic 
interaction with an environment, that, quite generally, 
can be thought of as a boson field. The precise form of the 
operator V which causes this behavior would, in general, 
depend on the detailed model of the environment. Here 
we choose the simplest possible operator, whose matrix 
elements are (2l| 



{e^\V\e,) = 




if i = 0, and < j < M 
otherwise, 



(5) 



where je^) is an eigenstate of the ground-state KS Hamil- 
tonian in the absence of the interaction with the bath, e^ 
is the corresponding eigenvalue, and the upper limit M 
is a given integer representing the number of states we 
keep in the simulation. (In the present case we have kept 
M = 15 states.) The parameter t^ gives the timescale 
over which the decay will occur, with larger values of id 
leading to longer decay times [22]. In the following we 
have chosen i^ = 1 fs. The operator V defined this way 
ensures that the stochastic Schrodinger equation (JH) is 
independent of the magnitude of r. 

Clearly, the above operator reduces the projection of a 
wave-function from the states {|ei), |e2) . . . |eA/-i)}, and 
increases the projection onto the ground state |eo). Phys- 
ically, it describes energy relaxation and dephasing. 

The stochastic Schrodinger equation dH) preserves 
the ensemble-averaged wave- function normalization [6|. 
However, the normalization is not necessarily satisfied 
for any particular realization of £{t). In order to reduce 
the number of dynamical calculations to perform the en- 
semble average, we have explicitly re- normalized |0) at 
every time step. As we will see below, with this approx- 
imation the decay into the ground state is evident even 
after a single realization of ' 



23| 



We begin by considering the behavior of an ensemble 
of He"^ ions interacting with the environment represented 
by the operator ([5]) . We prepare the system with all ions 
in the ensemble in the 2s state, denoted by |ei), and 
then let the electrons evolve in time according to equa- 
tion (|4]) [2J]. Panel (a) of Fig. [T] gives the projections 
Pi{t) = \{ei\(j>{t))\'^, as a function of time for one partic- 
ular realization of £{t). We see that the projection Po{t) 
onto the ground state approaches one as time evolves, 
while the projections onto other states approach zero, 
indicating energy relaxation into the ground state. In 
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FIG. 1: (Color online) Panel (a): Stochastic evolution of the 
projections Pi onto the unperturbed states e^, for He"*", as a 
function of time. The black, red, and blue curves correspond 
to Po, Pi, and P2, respectively. Note that the projections P3 
and Pi coincide with P2, since all three states have p sym- 
metry. All other projections are vanishingly small. Panel (b) 
Same as for panel (a), but averaged over five different runs, 
each with a different seed for the random number generator. 
Panel (c): Same as for panels (a) and (b), except that the 
dynamics were calculated using the Lindblad master equation 
^. Panel (d): Same as for panels (a)-(c), except the dynam- 
ics were calculated using the wave-packet collapse methodol- 
ogy of equation ([T]) . 



order to demonstrate that this behavior is not due to the 
particular choice of seed in our random number gener- 
ator, we also plot the projections Pi{t) averaged over 5 
different simulations with different seeds. One can clearly 
see that the fluctuations in panel (b) of Fig. [T] are reduced 
in comparison to the fluctuations in panel (a). 

For the single-electron case of He+ , we can analytically 



treat the ensemble average over all realizations of £{t) 
by considering the density matrix of this mixed state. 
Using equations ^ and ([3]), it can be shown [a, 0| that 
the resultant density matrix p evolves according to the 
Lindblad master equation 






rVpV'' -^pV^V-^V^'Vp. (6) 



In panel (c) of Fig. [T] we plot the matrix elements [ei 
showing the same behavior obtained with Eq. (|4]) [2£ 

We now discuss this result in terms of measurement 
theory. It is well-known that is possible to interpret the 
interaction with an environment as a continuous "mea- 
surement" of the state of the system - or, equivalently, 
of the state of the environment - with consequent non- 
unitary wave-packet reduction [8, 9]. We can make this 
point even clearer by assuming that every time the sys- 
tem interacts with the environment it emits a boson exci- 
tation (whether a photon or a phonon) and thus there is a 
finite probability dp — dt|(ei|(/))p/i(i that the the emitted 
excitation be detected by an apparatus [10| . Upon detec- 
tion of this excitation, the wave-function |0) collapses to 
the ground state |eo). This is the well-known postulate 
of wave- packet reduction. 

We can write the above in the form of a Schrodinger- 
type equation of motion that includes a stochastic vari- 
able 7(t), which has a probability distribution uniformly 
distributed between and 1. If 7 > dp, an emitted ex- 
citation is not detected, while if 7 < dp, the emitted 
excitation is detected, and the wavefunction collapses to 
the ground state. That is, during a small time At, |(/)) 
evolves according to 



m + At)) = e-'«-^*0(7(i) ~ dpMt)} 
+ 0(dp-7(i))|eo), 



(7) 



where ^(a;) is the Heaviside step function [2^. A similar 
approach has been used by Dalibard et al. in the context 
of quantum optics [10| . 

In Fig.[T]Jd), we plot the results from the time evolution 
of equation ([7]) for the problem of IIe+ relaxation, where 
we have considered only the ground state and the first 
excited state. Starting with \(j){t = 0)) = |ei), we evolved 
equation ([7|) in time for 1000 different realizations of 7(i), 
and found the average value of | (e^ | (/)(i) ) p , which we have 
denote by Pi. Note that each individual wave- function 
starts in the excited state, and then suddenly drops to 
the ground state the first time that 7(i) < dp. This wave- 
packet reduction occurs at a different time for each run, 
and the "remaining" excited states become exponentially 
less likely as time goes on; therefore, the average curve 
approaches an exponential. By comparing Fig.[TJd) with 
FigHJc) we see that the non- unitary wave-packet reduc- 
tion evolution is qualitatively similar to that provided 
by the stochastic Schrodinger equation ^, the differ- 
ence being in the fact that we have included only two 



states in the analysis of Eq. ([7]) . This equivalence there- 
fore illustrates a point of contact between the stochas- 
tic Schrodinger equation (j4]) and quantum measurement 
theory: the environment "measures" the state of the sys- 
tem, and, as a result, the wave-function is modified in a 
non-unitary way. 

We conclude by discussing the decay of neutral He, 
where a closed form for the KS single-particle Lindblad 
equation cannot be obtained. We prepare the system 
in such a way that both electrons (with spin a) are in 
the first excited state of the ground-state Hamiltonian, 
|(/)cr(i — 0)) = |ei). This means that the Pauli exclusion 
principle is automatically satisfied by our environment 
operator ([5]). For this case we compare the stochastic 
evolution with the one in which V = [27| . 

In Fig. ([2]), we plot the projections Pi = |(ei|0)p for 
the unitary evolution (panel (a)), as well as the projec- 
tions for the non-unitary evolution for one realization of 
£{t) (panel (b)) and averaged over five different realiza- 
tions of £{t) (panel (c)). As expected, in the presence 
of the environment, the projection onto the ground state 
I Co) approaches 1, while the occupations of other states 
are suppressed as time goes on. Here, however, we also 
note another effect of the interaction. Fig. [2ja) illustrates 
that, in the unitary evolution the projections Pi oscillate 
in time. This oscillatory behavior reflects the motion 
of the electrons as they alternately fall toward the nu- 
cleus, and then rebound outward [28|. Interaction with 
the environment has the effect of not only dampening 
these oscillations but also of modifying their frequency, 
the details of which vary depending on the particular re- 
alization of £{t). The introduction of the bath mediates 
new transitions for the single-particle wavefunction \(j)). 

Similarly to the case of He+, we can make a connec- 
tion with quantum measurement theory and study the 
decay of neutral He using Eq. ([7]). This is illustrated in 
Fig. [21(d). Here, however, we observe an important qual- 
itative difference. In the wave-packet reduction formal- 
ism described by Eq. ([7]), after the apparatus "detects" 
the excitation, the system immediately collapses onto the 
ground state |eo)) irrespective of the interactions of the 
other states with the environment. In contrast, the non- 
unitary evolution of the stochastic Schrodinger Eq. ^ 
involves a constant process of self- consistent interaction 
with the reservoir. This implies that the frequency of 
small oscillations is unchanged in the wave-packet reduc- 
tion formalism, while they change in time during dy- 
namical interaction with the environment as described 
by Eq. ©. Clearly, one could introduce this effect into 
the wave-packet reduction Eq. (O, but at a non-trivial 
complexity cost, while the stochastic Schrodinger Eq. ([6]) 
contains it naturally. 

In summary, we have used Stochastic TD-CDFT to de- 
scribe the interaction of an excited quantum system (He) 
with an external environment and its consequent decay 
into the ground state; a problem previously inaccessible 
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FIG. 2: (Color online) Panel (a): Projections Pi = |(ei|0>P 
for neutral He, as a function of time, for the case where the 
stochastic terms are not included (unitary evolution). The 
black, red, and blue curves correspond to Po, Pi, and Pg, 
respectively. Panel (b): Same for as panel (a), but with the 
inclusion of the interaction with the environment. Panel (c): 
Same for as panel (b), but averaged over five different runs, 
each with a different seed for the random number generator. 
Panel (d): Same as for panels (a)-(c), except the dynamics 
were calculated using the wave-packet collapse methodology 
of Eq. 0. 



via standard DFT methods. We have made a connec- 
tion of this open quantum problem with quantum mea- 
surement theory thus showing that Stochastic TD-CDFT 
may find applications in quantum information theory of 
realistic systems. 

We thank Y. Pershin for useful discussions and for 
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05ER46204. 
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In the calculations of this paper, we use small but finite 
time steps of duration At. The mean and autocorrelation 
of i{t) are thus given by I(ti) = and l{ti)l{tj) = Sijr/At, 
respectively, where ti and tj are arbitrary times, and 5ij 
is the Kronecker delta. We have chosen the probability 
distribution P{£) of £{t), to be Gaussian, so that P{£) = 

y/At/2TTTexp{-£^At/2T). 

In our numerical work we have used a position ba- 
sis so that the action of this operator in this ba- 

sis is {T\V^it)) = E,(rk.)E,(^.l^k.)fel<^W) = 

In Ref. 0, a similar operator for a two-state system is pro- 
vided to describe energy relaxation and decoherence. 
Nevertheless, this creates slight differences in the long-time 
limit of the dynamics with respect to the exact solution 
(see, e.g, Fig.[TJ. 

We represent the He nucleus with a simple 2/r poten- 
tial. We integrate out the singularity at the origin using 
a method similar to the Ewald method [ll|, [l^]. We use 
the Hockney method to calculate the potential of an iso- 
lated system ^^. The supercell is a cube of length 16.93 
A, and the grid spacing is 0.239 A. We use the Chebyshev 
method of constructing the propagator [14l.[l5J|. and we use 
a time step of 0.02 fs. 

Note that the density matrix procedure we have used here 
does not take into account the additional re-normalization 
that we apply in order to speed convergence in the cal- 
culations with Eq. Q. The close correspondence between 
panels (b) and (c) in Fig. [l] demonstrates that our forced 
normalization procedure does not cause significant devia- 



tion from the true average behavior of the system. 

It can be shown [10| that by defining an operator S such 

that S\ei) = |eo), the density matrix p associated with 

Eq. ([7]) satisfies the Lindblad equation ((Gj with S replacing 

V. 

We use the local density approximation to the scalar 

exchange-correlation potential [2, Isl] , as derived by Ceper- 



ley and Alder [l6l ] and parametrized by Perdew and 
Zunger [13]- For the current-density functional, we use the 
interpolation formula of Conti and Vignale [la |. 
Sugino et al. report similar high-frequency oscillations 
while studying the dipole moment of an isolated aluminum 
dimer [19(1. 



